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We study numerically the length distributions of the infrared monopole clusters 
" in pure SU(2) QCD. These distributions are Gaussian for all studied blocking steps 

of monopoles, lattice volumes and lattice coupling constant. We also investigate 
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the monopole action for the infrared monopole clusters. The knowledge of both the 
length distribution and the monopole action allows us to determine the effective 



entropy of the monopole currents. The entropy is a descending function of blocking 
scale, indicating that the effective degrees of freedom of the extended monopoles are 
getting smaller as the blocking scale increases. 



PACS numbers: 11.15.Ha,12.38.Gc,14.80.Hv 



INTRODUCTION 



The dual superconductor picture [l[ of the QCD vacuum is one of the most promising 
approaches to the problem of color confinement. This picture is based on the existence of 
Abelian monopoles in the vacuum of QCD. The monopoles are identified with the help of 
the so-called Abelian projection method 0, which is based on the partial gauge fixing of 



the SU(N) gauge symmetry up to an Abelian subgroup. The monopoles naturally appear 
in the Abelian projection due to compactness of the residual Abelian group. 

There are various numerical indications that the monopoles are responsible for the con- 
finement of quarks (for a review, see Ref. [3]). One of the most important observations is 
the monopole condensation in the low temperature (confinement) phase 0,0]. According to 
the dual superconductor mechanism the monopole condensation give rise to the formation 
of the chromoelectric string which confines the fundamental color sources. This expectation 
is confirmed by the fact that the non-zero tension of the chromoelectric string is dominated 
by the Abelian monopole contributions 0, 0, |^ . 

In the numerical simulations one observes that the trajectories of the Abelian monopoles 
form clusters, which can be divided by two ensembles: finite-sized clusters and one large 
percolating cluster 0, [1(1 11 1. The percolating cluster (or, infrared cluster) occupies the 



whole lattice while the sizes of the other clusters have an ultraviolet nature. The existence 
of the IR cluster is related to the monopole condensation This is understandable on an 
intuitive level, since generally the condensation is a microscopic effect reflecting itself in a 
zero-momentum component of the condensed field. The importance of the IR cluster for 
the confinement of quarks was also stressed in numerical calculations [10]: the tension of 
the confining string gets a dominant contribution from the monopoles belonging to the IR 
cluster, while the contribution of the UV clusters to the string tension is negligible. In the 
deconfinement phase the IR cluster disappears 0,113], as expected. 

In this paper we mostly concentrate on the numerical investigation of the properties of 
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the infrared monopole cluster. The length distributions and other properties of the UV and 
IR clusters were studied previously in Refs. In this publication we are 

going to investigate thoroughly the properties of the length distributions of the monopole 
clusters for various lattice volumes and sizes of the extended monopoles. 

The plan of the paper is the following. In Section |H] we describe the model and pro- 
vide the details of numerical simulations. Section IIHI is devoted to the investigation of the 
Abelian monopole action obtained by the inverse Monte-Carlo method. The distribution 
of the cluster length in the infrared clusters is studied in Section IIVI The knowledge of 
the monopole action and cluster distribution allows us, for the first time, to calculate the 
entropy of the lattice monopoles of various sizes. Our conclusions are presented in the last 
Section. 



II. MODEL AND SIMULATION DETAILS 

We study the pure SU(2) gluodynamics with the lattice Wilson action, S(U) = — |TrC/p, 
where (3 is the coupling constant and Up is the SU(2) plaquette constructed from the link 
fields. All our results are obtained in the Maximal Abelian (MA) gauge ^| which is defined 
by the maximization of the lattice functional 

R = J2^(^U(s,fi)a 3 U\s,fx)y (1) 

s,fi 

with respect to the gauge transformations U(s,fi) — ► £/(s,/i) = fl(s)U(s, (s + jl). The 
local condition of maximization can be written in the continuum limit as the differential 
equation {d ll + igA^){A l ^ — iA 2 ^) = 0. Both this condition and the functional ([!} are invariant 
under residual U(l) gauge transformations, f2 Abcl (u;) = di&g(e luj ( s \ e~ lu ^). 

The next step is Abelian projection of non-Abelian link variables to the Abelian ones 
after the gauge fixing is done. An Abelian gauge field is extracted from the 577(2) link 
variables as follows: 

~ f(l-\c(s,v)\ 2 ) 1/2 -c*(s,») \(u(s,fz) \ 

( '^"V (l-\c(s,ri\ 2 )V 2 ){ u*(8,ri)> [Z) 

where u(s,fi) = exp(i0(s, fi)) represents the Abelian link field and c(s, //) corresponds to 
charged matter fields. 

The Abelian field strength 6^ v (s) G (— 47r,47r) is defined on the lattice plaquettes by a 
link angle 9(s, fi) E [— ir, ir) as 9^ u (s) = 9(s, ji) + 9(s +/t, v) — 9(s + £>, fi) — 9(s, u). The field 
strength 9^(s) can be decomposed into two parts, 

9f JL y{s) = 9^(s) + 2imi^ u (s) , (3) 

where 9^(s) G [— n, n) is interpreted as the electromagnetic flux through the plaquette and 
m M ^(s) can be regarded as a number of the Dirac strings piercing the plaquette. 

The elementary monopole currents is conventionally constructed using the DeGrand- 
Toussaint 3 definition: 



k v( s ) = ^ e ^ P ad u m pa (s + fi), (4) 
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where d is the forward lattice derivative. The monopole current is defined on a link of 
the dual lattice and takes values 0, ±1,±2. Moreover the monopole current satisfies the 
conservation law automatically, 

&Ms) = 0, (5) 

where d 1 is the backward derivative on the dual lattice. 

Besides the elementary monopoles one can also define the so called extended 
monopoles 0. In this paper we use the type-2 construction according to the classifica- 
tion of the extended monopoles adopted in Ref. 0. The n 3 extended monopole is defined 
on a sublattice with the lattice spacing b = na, where a is the spacing of the original lattice. 
Thus the construction of the extended monopoles corresponds to a block spin transformation 
of the monopole currents with the scale factor n, 

n-l 

k^\s)= k^ns + (n-l)ii + ii> + jp + la) . (6) 

i,j,l=0 

The Abelian dominance and the monopole dominance in the infrared region of QCD 
implies that at least important infrared observables (such as the fundamental string tension) 
can be calculated using the Abelian fields or the monopole degrees of freedom only. 

In what follows we discuss an effective model of the monopole currents corresponding to 
pure SU(2) QCD. Formally, we get this effective model through the gauge fixing procedure 
applied to the original model. Then we integrate out of all degrees of freedom but the 
monopole ones. An effective Abelian action is related to the original non-Abelian action 
S[C, 9} (matter, C, and Abelian gauge, 9, fields, Eq. (J2J)) as follows: 

Z = JvuU DCe~ s ^S(X) A FP (U)j = J Vu e~ s ^ . (7) 

Here and below we omit irrelevant constant terms in front of the partition functions. In 
Eq. (J7J) the term 8(X) represents the gauge-fixing condition 1 and A FP (U) is the correspond- 
ing Faddeev- Popov determinant. Next step is to relate the effective monopole action to the 
effective U(l) action: 



oo „ 

n £ ) V98(k,(s)-k,(s;9)))e- s ^m (8) 

S >M fc M 0)=-oo 
oo 

n e )(nwO e ~ larw ' < 9 > 

Av(s)=-oo s 

where k^(s;9) is the monopole current defined as a function of the Abelian fields, 9, via 
relations (j2J) and (J3J). 

Our simulation statistics is represented in Table The gauge configurations were gen- 
erated with the help of the standard Monte-Carlo algorithm. In most simulations we use 



1 As we have discussed above, the MA gauge fixing condition is given by a maximization of the functional 
and therefore the use of the local condition X = 0, implied in Eq. JJJ), is a formal simplified notation. 
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Lattice 


P 


Blocking 


Configuration 


size 




factor 


number 


6 


2.1~2.4 


1 


3000 


8 


2.1~2.4 


1 


3000 


10 


2.1~2.4 


1 


3000 


12 


2.1~2.4 


2 


3000 


1 A 

14 


2.1~2.4 


1 


3000 


16 


2.1~2.4 


2 


3000 


24 


2.1~2.4 


2,3,4 


3000 


32(SA) 


2.1~2.6 


2,3 


950 


48 


2.1~2.6 


2,3,4,6,8 


2200 



TABLE I: Simulation statistics. 

the usual iterative algorithm to fix the MA gauge. However, in order to check the Gribov 
copy dependence of the MA gauge fixing we also use the so called simulated annealing (SA) 
algorithm with five Gribov copies. We refer a reader for a detailed description of the SA 
method to Ref. [T3|, where the advantages of the SA method compared to the iterative 
algorithm are illustrated. 



III. MONOPOLE ACTION FOR VARIOUS CLUSTERS 



It is well known that in gluodynamics the monopole trajectories can be separated into 
the infrared and ultraviolet monopole clusters. There is only one IR monopole cluster which 
occupies all volume of the lattice, and a large number of shorter monopole trajectories (UV 
clusters). In Figure ^ we show the typical length distributions, D(L), of the monopole 
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FIG. 1: Distribution of the lengths of the monopole trajectories at various [3 for (a) elementary 
and (b) n = 2 blocked monopoles. 
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trajectories in all clusters. We show the data for elementary, n = 1, and blocked, n = 2, 
monopoles at various lattice coupling constants (3. One can see that for all considered values 
of the coupling (3 the infrared cluster and the ultraviolet clusters can be unambiguously 
separated due to a wide gap between them. Moreover, the distributions of the elementary 
and blocked monopoles are qualitatively similar. 

Note that at zero temperature the gap between IR and UV clusters becomes smaller as 
the physical lattice size decreases. This behaviour can be observed in Figure^ At very small 
lattice size the gap between UV and IR clusters disappears and the IR and UV clusters can 
not be resolved. This behaviour of the monopole clusters leads to the deconfming transition 
("crossover") which takes place in sufficiently small physical volumes. 

The distribution of the ultraviolet clusters was studied both numerically (ill l3 | and 



analytically |13j, LLSj]- The distribution can be described by a power law D uv oc L r , where 



| . This behaviour indicates that the monopoles in UV 
3. In our simulations we are mainly concentrated on 



the power r is very close to 3, Ref. 
clusters show a random walk picture 
the IR monopole cluster because, as we have already mentioned in the Introduction, the IR 
cluster is important for the confinement of quarks. Below we study the monopole-related 
quantities using the largest monopole cluster only unless stated otherwise. 

In general, the monopole action, S™ n , can be represented as a sum of the n-point (n > 2) 
operators Si, Ref. [3. 

S[k] = Y^fiSi[k], (10) 

i 

where /j are coupling constants. In this paper we adopt only the two-point interactions 
in the monopole action {i.e. interactions of the form Si ~ k tl (s)k IJ ,'(s > )). Following Ref. [H 
we derive the effective monopole action (fTU|) from the configurations of monopole currents, 
{fc M (s)} using an inverse Monte-Carlo method. The original monopole configurations were 
generated by the usual Monte-Carlo algorithm of SU(2) gluodynamics. 

The dominant term in the monopole action (jl(J|) corresponds to the most local self- 
interaction of the monopole currents, Si[k] = 5^ s „A^(s). The contributions to the action 
associated with other interactions are small compared to the leading term. As an example 
we show the leading contribution and the full action associated with the IR monopole cluster 
for (3 = 2.4 and n — 1, 2 in Figure 121 Moreover, one can find that both the monopole action 
and the self-coupling contribution to it are proportional with a good accuracy to the length 
of the monopole loop. 

In Figure E3 we plot the ratio S[k]/L for various lattice volumes and blocking sizes 2 . One 
can notice that the coupling constant depends on the product b = a ■ n and almost does 
not depend on the variables a and n separately, in agreement with observations of Ref. ■ 
Below we will observe this type of scaling in many other monopole quantities. Another 
observation is that the monopoles obtained with the SA procedure have the same action as 
the monopoles defined in the MA gauge which is fixed by the usual iterative algorithm. 

It would also be interesting to compare the monopole action associated with the IR 
cluster and the action associated with the whole monopole ensemble. The simplest quantity 
to compare is the fx self-coupling parameter which is a dominant coupling in the action. 
In Figure E] we show j\ for both ensembles. First, we easily notice that for chosen lattices 



2 In this figure and all other figures below we plot all dimensional quantities in units of the string tension, a . 
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FIG. 2: The total monopole action and the contribution of the self-interaction term to the action 
for (a) elementary and (b) n = 2 extended monopoles vs. length of the monopole trajectory in the 
IR cluster as calculated on 24 4 lattice at (3 = 2.4. 
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FIG. 3: The ratio S/L, in physical units, as the function of b for various lattices, N 4 , and blocking 
steps, n. 



the coupling constant fi is independent of the lattice volume. Second, we see that for large 
blocking scales b the type of the ensemble (the IR cluster or the whole ensemble) is not 
essential for determination of f\. However, at small b values, by/a < 0.5, the type of the 
lattice ensemble becomes important, since in this region 

f m > jtoui j for ^< .5. (11) 

The observed difference between the couplings can be affected by finite-size effects since the 
leftmost points in our data correspond to elementary (of size a) monopoles. Moreover, in our 
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FIG. 4: The self-interaction coupling constant fi as the function of b calculated for the largest 
monopole cluster and for the whole monopole ensemble on lattices 24 4 and 48 4 . 



studies we have included only the two-point interactions in the monopole action (jlOjh while 
the two-point action becomes unreliable at too small values of b, and one has to include 
higher-point interactions. 

Despite a possible influence of the lattice artifacts, the observation pip may have a phys- 
ical meaning related to the simple fact that the larger coupling f\ the smaller density of the 
monopoles is. Thus Eq. (fTTj) is in agreement with the numerical fact at large lattice cou- 
pling (3 (i.e., at small lattice spacing a) the density of the monopoles in the largest cluster 
is noticeably smaller than the total monopole density Figure 0] is also in a qualita- 

tive agreement with the fact 2l| that the excess of the Abelian action around elementary 
monopoles in all clusters and in the IR cluster almost coincide with each other. However, 
the larger the physical size of the monopole cube the better the agreement between the 
all-cluster and IR-cluster actions is, in accordance with Figure EJ 



IV. MONOPOLE LENGTH DISTRIBUTION FOR IR CLUSTER 

As we have mentioned, the length distribution of the monopole trajectories in the ultra- 
violet clusters was found to obey the power-law. It would also be interesting to study the 
length distribution for the infrared clusters following Ref. ^(J . 

Since the density of the elementary monopoles from infrared clusters is finite (in terms of 
physical units) in the continuum limit [20|, we may expect that the density of the extended 
monopoles (with a fixed blocking scale b) is finite as well. The finiteness of the density is 
consistent with the observation that the monopole length distribution is localized around a 
certain value of the monopole length, L max (see Figure P). This value should be proportional 
to the physical volume, V, of the system, L max cx V. Indeed, as one can qualitatively judge 
from Figure ^ the position of the peak of the IR length distribution increases with increase 
of the physical volume of the system (i.e., with decrease of the lattice coupling (3). 
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The length distribution function, D(L), is proportional to the weight with which the 
particular trajectory of the length L contributes to the partition function. On the other 
hand, the action of a monopole trajectory is proportional to the length of the trajectory, 
S oc L, as we have illustrated in the previous Section. Thus the monopole action contributes 
in a form of an exponential factor, oc e~f L , to the weight with which this trajectory appears in 
the partition function. Here / is a parameter which is close to the self-coupling fi according 
to FigureEJ The entropy of the monopole trajectory also contributes to the monopole length 
distribution, which is proportional to \i L (with fi being positive number) for sufficiently large 
monopole lengths, L. Thus the distribution of the monopole trajectories in infinite volume 
must be described by the function 

D!J}(L)<xn L -e-f L = er L , -y = hxfi-f. (12) 

In this equation we neglect a power-law prefactor which is essential for the distribution of 
the ultraviolet clusters 3 3 • 



The observed localization of the infrared cluster distribution imply that the finite volume 
provide a certain cut which depends on the volume of the system. The simplest distribution 
of this kind can be described by the function: 

D IR {L) = exp{-aL v + 7 L} , (13) 

where a, 7 and rj are certain parameters. As we find below, the parameter rj which character- 
izes the cut due to the volume effect, is close to 2 with a big accuracy, 77 ~ 2. Moreover, as we 
mentioned, the parameter 7 is characterizing the action and the entropy of the monopole cur- 
rents and thus it should depend only of the physical size of the blocked monopole, 7 = 7(6). 
As for the parameter a, it should also be dependent of the volume of the system, a = a(b, V). 
Thus we employ the following parameterization of the IR monopole distribution at finite 
volume: 

D IR {L) = exp{-a(6, V)L 2 + 7 (6)L} . (14) 
The peak of the distribution (JHj), 

7(6) 



2a(6, V) 



(15) 



is expected to be proportional to the volume of the system, L max oc V to insure the finiteness 
of the IR monopole density, 

L m ax _ 7(^) (-\a\ 

PIR - — - 2a(b,V)V (16) 
in the thermodynamic limit, V — > 00. Thus from Eq. (jl6|) we conclude that 

a(b 1 V) = A(b)/V } (17) 



Below we work with the distribution of the pure exponential form H12(l . We also repeated our analysis 
with the prefactor L~ 3 included. We observed that the results with and without the power-law prefactor 
are the same within the small error bars. 
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FIG. 5: The original histograms of the length distribution in the IR cluster are shown by grey color. 
The averaged distributions are shown by circles, and the fits by the function (|14j) are represented 
by the dashed line. 



where the function A(b) depends only on the size of the blocked monopole, b. Eq. (jl7j) 
implies that in the thermodynamic limit the parameter a vanishes and the finite-volume 
distribution (JT3J) is reduced to Eq.(JT2J), as expected. 

We show typical examples of the IR cluster distributions in Figure 03 by grey color. One 
can see that these histograms have an almost symmetric structure, but due to the lack of 
statistics these histograms can not be fitted by the function (I14j) . In order to show that the 
distribution of the monopole follows Eq. (jl4|) we smooth the data by increasing the step of 
the histograms (which was equal to 2) and then averaging the data inside the coarse steps. 
We show the averaged (and suitably rescaled) histograms and their fits by the function (JT3J) 
in the same Figure. One can see that the averaged histograms are very close to the Gaussian 
distribution. Similar behaviour can also be observed for all IR monopole cluster distributions 
we have studied in this paper. 

In order to justify the chosen value of the parameter 7] in Eq. (fT^) we have also fitted the 
averaged histogram data by Eq. (|T3~|) in which rj is treated as a fitting parameter. The best 
fit (for (3 = 2.4 and n = 1 on 24 4 lattice, as an example) gives us the result r] = 2.05(15). 
Fits of other histograms give us similar results. Thus we fix below r] = 2. 

The histograms in Figure were obtained with a high simulation statistics (3000 con- 
figurations according to Table HJ). In order to get a perfect gaussian we would need much 
more statistics which would consume a lot of CPU time. To avoid this we assume that the 
numerical data for length distribution of the IR monopoles is described by Eq. (jl4|) . This 
assumption is justified by the analysis we have performed above. Then one can evaluate the 
central values of the parameters a and 7 using the simple formulae (valid for a Gaussian 
distribution): 

_ 1 1 (L) 
2 <L 2 ) - (L) 2 ' 7 (L 2 )-(L) 2 ' 

where the averaging (■ • •) is performed using weights from the histograms. 

To evaluate the errors for the parameters a and 7 we use the standard bootstrap method. 
Namely, we make a resampling of the original data describing lengthes of the IR monopole 



(18) 
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clusters, L max . We construct a resampled configuration by selecting n con f random values of 
L max (note that a single value of L max can be picked up multiple number of times), where 
n con f is the total number of the monopole configurations. Then we evaluate the values of 
a and 7 at each resampled configuration using Eq. (JTHJ). The distribution of these values is 
again the Gaussian with the width equal to the corresponding error. We plot examples of 
the histograms for a and 7 values in Figures |Bfa,b). 




(a) (b) 

FIG. 6: The distribution of the parameter a and 7 for elementary monopoles at j3 = 2.4 on 24 4 
lattice. The fits by a gaussian function are shown by the solid lines and the value of the errors are 
indicated by shadowed regions. 

We have checked the applicability of Eqs. (fTHJ) and the use of the bootstrap method 
on a smaller, 16 4 , lattice. Namely, we have generated length distributions using from low 
statistics (2000 configurations) to high statistics (10 5 configurations) ensembles. We used 
the bootstrap method along with Eqs. (fTHjl to evaluate the coefficients a and (3 for the 
distribution measured with the lowest statistics. On the other hand, the high statistics 
distribution is a (almost perfect) Gaussian and therefore we get the desired coefficients 
directly from the fit (JHJ). The comparison of the coefficients shows that the central values 
as well as the estimated errors for the low and for the high statistics ensembles coincide with 
each other within a few percents. We illustrate our analysis in Figure [3 for /? = 2.1 and 
(3 = 2.2 using the parameter 7 as an example. The values of 7 obtained with the standard 
method are plotted vs. number of configurations, N con f, used in the analysis. The horizontal 
lines represent the results coming from the bootstrap method applied to the low-statistic 
ensemble (the statistical errors are indicated by shadowed regions). We conclude that the 
bootstrap method allows to get reliable results using the distributions with low-statistics. 

In order to confirm our expectation (|17|) we plot the parameter a vs. the ratio N/n in 
Figure |Sfa) for selected set of coupling constants (5 and the blocking steps of the monopole, 
n. Since the volume of the blocked lattice is (N/n) 4 , we expect that the parameter a behaves 
as a oc (N/n)~ 4 . This behaviour is seen in Figures IHJa,b). The parameter a multiplied by 
the lattice volume almost does not depend on the lattice size N according to Figure E[b). 

According to our discussion above the fitting parameter 7 should only be a function of 
the blocking size b and should not depend on the volume of the lattice. In Figure^ we show 
the parameter 7 is indeed independent of the lattice size N. 
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FIG. 7: Check of the bootstrap method on 1643 lattice using the parameter 7 as an example (the 
explanation is in the text). 
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FIG. 8: (a) The fitting parameter a as the function of the size N/n of the coarse lattice; (b) the 
parameter a multiplied by the lattice volume as the function of the lattice size N. 



The fitting parameters a and 7 are shown as functions of the physical scale b in Fig- 
ures a) and (b), respectively. The parameter 7 shows the scaling behaviour in a sense 
that it depends on the blocking step n and lattice spacing a in the form of the product 
b = n ■ a. 
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FIG. 9: The illustration of the independence of the fitting parameter 7 on the lattice size, N. 
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FIG. 10: The fitting parameters (a) a and (b) 7 as functions b for various lattice volumes, iV 4 , and 
monopole blocking steps, n. 

V. MONOPOLE DENSITY AND ENTROPY 

A. Monopole density 

The simplest physical characteristic of the monopole ensemble is its density. It is inter- 
esting to compare the monopole density obtained from the IR monopole cluster distribution, 
Eq. (|lb|). with the direct observation of the monopole density, 



Pir 



4(na) 6 ■ (N/n) 



4 \ ' -> 1 M 



\k { :\s)\). 



(19) 
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Here the blocked monopole current, kj?\ is defined by Eq. (JBJ). The normalization factor in 
Eq. ()19j) appears naturally if one notes that b = na and 4(iV/n) 4 are the lattice spacing and 
the number of links of the coarse lattice, respectively. 

If the fitting function ()14J) describes the data correctly then one should observe no essential 
difference between the infrared monopole density obtained from the fits of the monopole 
distributions, (|14I16|) compared to the density obtained in a direct way (fK?|) . This is indeed 
the case according to Figure HIT a). 

Another information which can be extracted from this Figure is that the blocked monopole 
density goes to a fixed limit, lim^o P ~ 0.9cx 3//2 , as the blocking size b gets smaller. It is 
also possible that the monopole density shows a wide plateau around by/a ~ 0.2. In order 
to discriminate between these options one should study the blocked monopole density at 
smaller values of lattice spacing, a, and, consequently, at larger lattice volumes. We also 
note that the value of the blocked monopole density quoted above is about 30% larger than 
the value of density 2(| of the elementary infrared monopoles in the continuum limit. 

The monopole density is known to be sensitive to the details of the gauge fixing pro- 
cedure ■ In order to check the effect of the gauge fixing we compare in Figure ITTT b) 
the infrared monopole density obtained using the SA and iterative gauge fixing algorithms. 
One can see from this Figure that at large b there is practically no difference between the 
monopole densities obtained with the use of the different algorithms. However, there exists 
some difference at small b since the SA monopole density is smaller than the density ob- 
tained with the help of the iterative algorithm. This slight dependence of the density on the 
gauge fixing algorithm at small b may explain the discrepancy between our results and the 
results of Ref. [2(j] mentioned above. Another source of the discrepancy is the qualitative 
difference between the elementary and the blocked monopoles. Since the scale b is taken to 
be independent of the lattice spacing a while a tends to zero in the continuum limit, the 
elementary monopoles are expected to be more affected by the ultraviolet lattice artifacts. 
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B. Monopole entropy 

The distribution of the monopole trajectories depends both on the monopole action and 
on the monopole entropy as we have already discussed in Section IIVI Therefore the knowl- 
edge of the distribution and the monopole action allows us to define the entropy of the 
monopole currents. If the monopoles make a simple random walk on the four-dimensional 
hypercubic lattice then the entropy factor for elementary monopoles is expected to be equal 
to seven, /i = 7, since there are seven choices at each site for the monopole current to go 
further (the monopole trajectory is obviously non-backtracking due to the presence of the 
magnetic charge). 

The balance between energy and entropy of the elementary monopole trajectories plays 
an important role. For example, the compact U(l) gauge model in four dimensions possesses 
a phase transition associated with the monopole (de) condensation which is defined as a point 
on the phase diagram where the entropy and the energy of the monopole trajectories are 
the same. This point is located by the condition 7 = 0, where 7 characterizes the monopole 
distribution, Eq. (fT5|) or (ffl jl . 

In the case of the compact U(l) gauge model the coefficient / in the action of the 
monopole trajectory, S = fL, is proportional to the lattice gauge coupling, fu^ oc Pu(i) and 
fi = 7. This fact allowed authors of Ref. |22| to find the critical value of deconfinement 
transition point with a great accuracy. Indeed, if 7 is negative then the infrared cluster 
disappears and the confinement of electric charges is lost. The energy-entropy balance was 
also studied numerically for the monopoles in compact U(l) gauge theory |23| and in finite- 
temperature pure SU(2) gauge theory fib] ]. 

In the pure zero temperature QCD the coupling 7 is positive at all values of the lattice 
coupling constant 4 (3. The approximate cancellation of the entropy factor and the energy 
of the elementary monopoles in the zero-temperature gluodynamics is known as "fine tun- 
ing" k| 3|. This fact is in agreement with the existence of physical scaling of the infrared 
cluster in zero-temperature pure QCD. 

The entropy factor /i of the infrared monopole trajectories can be obtained from the IR 
cluster distribution and the monopole action according to Eq. (fT2"|) . 

fi = e 7+/ . (20) 

We calculate numerically the parameters 7 and / to find the entropy factor /i for various 
scales b and lattice sizes. Our results are presented in Figure H21 The entropy shows an 
approximate scaling behaviour in a sense that the entropy depends only on the scale b and 
is independent of the lattice spacing, a, and the blocking factor, n, separately. One can also 
notice that the entropy /i is independent of the volume of the lattice. The largest scaling 
violations happens at small blocking sizes n = 1,2 at which the finite-size artifacts are 
expected to be strong. 

The entropy factor fi is a declining function of the scale b. According to discussion 
above one can expect that for elementary monopoles the factor /1 should be equal to seven. 
However, /1 > 7 for small values of b, as can be seen from Figure We explain this 
small-6 behaviour as an artifact of our numerical procedure adopted in this paper. Indeed, 



4 The coupling (5 must be smaller then a certain value at which the unphysical deconfinement transition 
happens at finite volume. 
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FIG. 12: Entropy factor [i vs. b. The dashed line represents the fit by Eq. with free (/-parameter 
and the solid line corresponds to the fixed parameter, g = 1. 



we have used the quadratic monopole action while at small b higher-point interaction terms 
are essential and thus the monopole action can not be reliably described by the quadratic 
terms only [l9| . 

At large b the entropy factor (|20|) is smaller than seven. Formally this means that the 
motion of the blocked monopoles is constrained. We have fitted the entropy by the function: 

/^/ico + CV", (21) 

where /too, C and q are the fitting parameters. The best fit is shown in Figure H~2l bv the 
dashed line. The corresponding best fit parameters are: /i^ = 1.6(4), C = 1.7(5) and 
q = 1.2(2). The most interesting fitting parameter is /too which is the asymptotic value 
of the entropy in the infrared limit 6a 1 / 2 — > oo according to Eq. (J21|) . Unfortunately, the 
value of the asymptotic entropy is obtained with a big error bar in the above fit. In order 
to increase the accuracy we notice that the power q is very close to unity. Fixing /i = 1 in 
Eq. (j2Tj) and repeating the fitting procedure again we get /too = 1-15(25) and C = 2.2(1). 
The corresponding best fit curve is shown in Figure by the solid line. 

The fact that the asymptotic value of the entropy is very close to the unity in the limit 
bo 1 ' 2 — ► oo may have a simple explanation. The monopole with a large blocking size 
b behaves as a classical object and its motion is no more a simple random walk. The 
predominant motion of the large-6 monopole is close to a straight line. 



VI. CONCLUSION 

We studied numerically the distributions of the infrared monopole currents of various 
blocking sizes, n, on the lattices with different spacings, a, and volumes, N. The distributions 
can be described by a gaussian anzatz with a good accuracy. The anzatz contains two 
important terms: (i) the linear term, which possesses information about the energy and 
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entropy of the monopole currents; and (ii) the quadratic term, which suppresses too large 
infrared clusters. The linear term is independent of the lattice volume while the quadratic 
term is inversely proportional to the volume. The monopole density determined from the 
parameters of the gaussian fits coincides with the result of the direct numerical calculation. 

We also studied the action of the monopoles belonging to the infrared clusters and com- 
pared it with the action of the total monopole ensemble. It turns out that the self-coupling 
coefficients for both these ensembles are almost the same at large b. However, as the blocking 
scale b is decreased the self-coupling coefficient for infrared monopole cluster gets noticeably 
larger then the coefficient for the total monopole ensemble. This can be explained by the 
fact that the self-interaction coefficient depends on the monopole density (the larger density 
the smaller the coefficient is), and the difference between the total density and the infrared 
density increases as b gets smaller. 

The knowledge of both the coefficient in front of the linear term of the gaussian distribu- 
tion and the monopole action for infrared clusters allows us to determine the entropy factor 
of the extended (blocked) monopole currents. We have numerically shown that the entropy 
of the blocked monopole currents is a descending function of b = na, indicating that the 
effective degrees of freedom of the blocked monopoles are getting smaller as the blocking 
scale b increases. This corresponds to the classical picture: the monopole with the large 
blocking size b becomes a macroscopic object and the motion of such a monopole is close to 
a straight line. 
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